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We have studied the spin ordering of a dilute classical Heisenberg model with spin concentration x, 
and with ferromagnetic nearest- neighbor interaction Ji and antiferromagnetic next-nearest-neighbor 
interaction J2. Magnetic phases at absolute zero temperature T — are determined examining the 
stiffness of the ground state, and those at finite temperatures T 7^ are determined calculating 
the Binder parameter ql and the spin correlation length ^l- Three ordered phases appear in the 
x — T phase diagram: (i) the ferromagnetic (FM) phase; (ii) the spin glass (SG) phase; and (iii) the 
mixed (M) phase of the FM and the SG. Near below the ferromagnetic threshold a^p, a reentrant SG 
transition occurs. That is, as the temperature is decreased from a high temperature, the FM phase, 
the M phase and the SG phase appear successively. The magnetization which grows in the FM 
phase disappears in the SG phase. The SG phase is suggested to be characterized by ferromagnetic 
clusters. We conclude, hence, that this model could reproduce experimental phase diagrams of 
dilute ferromagnets Fea;Aui_a; and Eua:Sri_a;S. 

PACS numbers: 75.10.Nr,75.10.Hk,75.10.-b 



I. INTRODUCTION 



Prototypes of spin glass (SG) are ferromagnetic di- 
lute alloys such as Fea;Aui_2:0|, Eua;Sri_:rS|3, 1^ and 
Fea;Ali_a;P, Those alloys have a common phase di- 
agram as schematically shown in Fig. 1. It is shared 
with the ferromagnetic (FM) phase at higher spin con- 
centrations and the SG phase at lower spin concentra- 
tions, together with the paramagnetic (PM) phase at 
high temperatures. A notable point is that a reentrant 
spin glass (RSG) transition occurs at the phase bound- 
ary between the FM phase and the SG phase. That is, 
as the temperature is decreased from a high tempera- 
ture, the magnetization that grows in the FM phase van- 
ishes at that phase boundary. The SG phase realized 
at lower temperatures is characterized by ferromagnetic 
clusters jl|,i2i,|3i,@|- A similar phase diagram has also been 
reported for amorphous alloys (T'i_2:T^)75B6Al3 with T 
= Fe or Co and T' = Mn or Ni'tf]. It is believed that the 
phase diagram of Fig. 1 arises from the competition be- 
tween ferromagnetic and antiferromagnetic interactions. 
For example, in Fe^jAui-a;, the spins are coupled via the 
long-range oscillatory Ruderman-Kittel-Kasuya-Yoshida 
(RKKY) interaction. Also, in Eua;Sri_2;S, the Heisenberg 
spins of 5* = 7/2 are coupled via short-range ferromag- 
netic nearest-neighbor exchange interaction and antifer- 
romagnetic next-nearest-neighbor interaction'^. Never- 
theless, the phase diagrams of the dilute alloys have not 
yet been understood theoretically. Several models have 



been proposed for explaining the RSG transition[^l^[lo|. 
However, no realistic model has been revealed that re- 
produces it[ll[ll[3. Our primary question is, then, 
whether the experimental phase diagrams with the RSG 
transition are reproducible using a simple dilute model 
with competing ferromagnetic and antiferromagnetic in- 
teractions. 

This study elucidates a dilute Heisenberg model with 
competing short-range ferromagnetic nearest-neighbor 
exchange interaction Ji and antiferromagnetic next- 
nearest- neighbor interaction J2. This model was ex- 
amined nearly 30 years ago using a computer simula- 
tion technique 01 at rather high spin concentrations and 
the phase boundary between the PM phase and the FM 
phase was obtained. However, the SG transition and the 
RSG transition have not yet been examined. Recent ex- 
plosive advances in computer power have enabled us to 
perform larger scale computer simulations. Using them, 
we reexamine the spin ordering of the model for both 
T — Q and T ^ in a wide-spin concentration range. 
Results indicate that the model reproduces qualitatively 
the experimental phase diagrams. In particular, we show 
that the model reproduces the RSG transition. A brief 
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FIG. 1: A schematic phase diagram of a ferromagnetic dilute 
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report of this result was given in Ref. 15. 

The paper is organized as follows. In Sec. II, we 
present the model. In Sec. Ill, the ground state prop- 
erties are discussed. We will determine threshold a;F, 
above which the ground state magnetization remains fi- 
nite. Then wc examine the stabilities of the FM phase 
and the SG phase calculating excess energies that are ob- 
tained by twisting the ground state spin structure. Sec- 
tion IV presents Monte Carlo simulation results. We will 
give both the phase boundaries between the PM phase 
and the FM phase and between the PM phase and the 
SG phase. Immediately below x — xf, we find the RSG 
transition. Section V is devoted to our presentation of 
important conclusions. 
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FIG. 2: Ground state magnetizations AI in Lx Lx L lattices 
for various spin concentrations x. 



II. MODEL 

We start with a dilute Heisenberg model with compet- 
ing nearest-neighbor and next-nearest-neighbor exchange 
interactions described by the Hamiltonian: 



- JiXjXjSj ■ Sj + y^^ J2XkXiSk ■ S, 

(ij) (kl) 



(1) 



where Si is the classical Heisenberg spin of = 1; 
Ji(> 0) and J2(> 0) respectively represent the nearest- 
neighbor and the next-nearest-neighbor exchange inter- 
actions; and Xi — 1 and when the lattice site i is occu- 
pied respectively by a magnetic and non-magnetic atom. 
The average number of x{= (xi)) is the concentration of 
a magnetic atom. Note that an experimental realization 
of this model is Eu2:Sri_a;S[3|, in which magnetic atoms 
(Eu) are located on the fee lattice sites. Here, for sim- 
plicity, we consider the model on a simple cubic lattice 
with Ja = 0.2JiO. 



III. MAGNETIC PHASE AT T = 

We consider the magnetic phase at T = 0. Our strat- 
egy is as follows. First we consider the ground state of 
the model on finite lattices for various spin concentrations 
X. Examining the size dependence of magnetization M, 
we determine the spin concentration xp above which the 
magnetization will take a finite, non-vanishing value for 
L — > CX3. Then we examine the stability of the ground 
state by calculating twisting energies. We apply a hybrid 
genetic algorithm (HGA) [13 for searching for the ground 
state. 



A. Magnetization M at T = 

We treat lattices oi L x L x L with periodic bound- 
ary conditions. The ground state magnetizations Ai'^(= 
^^"6 calculated for individual samples and av- 
eraged over the samples. That is, M = [\M2\], where 
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FIG. 3: Binder parameter ql at T — in L x L x L lattices 
for various spin concentrations x. 



[■ ■ •] represents a sample average. Numbers of sam- 
ples with different spin distributions are Ng = 1000 for 
L < 8, Ns = 500 for 10 < L < 14, and A^^ = 64 for 
L > 16. We apply the HGA with the number Np of par- 
ents of Np = 16 for L < 8, A^p = 64 for L = 10, Np = 128 
for L 12, . . . , and Np = 512 for L > 16. 

Figure 2 portrays plots of magnetization M as a func- 
tion of L for various spin concentrations x. A consid- 
erable difference is apparent in the L-dependence of M 
between x < 0.82 and x > 0.84. For x < 0.82, as L in- 
creases, M decreases exponentially revealing that M ^ 
for L ^ 00. On the other hand, for x > 0.84, M de- 
creases rather slowly, suggesting that M remains finite 
for L "> 00. 

To examine the above suggestion, we calculate the 
Binder parameter defined as 



<7L = (5 - 3 



[|MG|2]2 



)/2. 



(2) 



When the sample dependence of M2 vanishes for L 
00, ql increases with L and becomes unity. That is, if the 
system has its magnetization inherent in the system, gL 
increases with L. On the other hand, (7l — > for L —>■ 00 
when M2 tends to scatter according to a Gaussian dis- 
tribution. Figure 3 represents the L-dependence of gL 
for various x. For x < 0.82, as L increases, increases 
and subsequently becomes maximum at L ~ 8, decreas- 
ing thereafter. This fact reveals that the FM phase is 
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absent for x < 0.82. For x > 0.84, a decrease is not 
apparent. In particular, for x > 0.86 increases grad- 
ually toward 1, indicating that the FM phase occurs for 
L — > c». We suggest, hence, the threshold of the FM 
phase of xp = 0.84 ± 0.02 at T = 0. 

B. Stiffness of the ground state 

The next question is, for x > xp, whether or not the 
FM phase is stable against a weak perturbation. Also, for 
X < xp, whether or not some frozen spin structure occurs. 
To consider theseproblems, we examine the stiffness of 
the ground state plll20l|. 

We briefly present the methodfigf. We consider the 
system on a cubic lattice with Lx Lx{L + l) lattice sites 
in which the z-direction is chosen as one for {L + 1) lattice 
sites. That is, the lattice is composed of {L + 1) layers 
with LxL lattice sites. Periodic boundary conditions are 
applied for every layer and an open boundary condition 
to the z-direction. Therefore, the lattice has two oppo- 
site surfaces: fii and f^L+i. We call this system as the 
reference system. First, we determine the ground state 
of the reference system. We denote the ground state spin 
configuration on the Ith layer as {<S^z,i} {I = 1 - {L + 1)) 
and the ground state energy as E2- Then we add a dis- 
tortion inside the system in such a manner that, under a 
condition that {Si^i} are fixed, {S'l+i.i} are rotated by 
the same angle around some common axis. We also 
call this system a twisted system. The minimum energy 
El{4>) of the twisted system is always higher than E'^. 
The excess energy Ai?L ((/))(= El{4>) — E^) is the net 
energy that is added inside the lattice by this twist, be- 
cause the surface energies of fii and Hl+i are conserved. 
The stiffness exponent 9 may be defined by the relation 
AEL{(j)) oc L^im. If 6* > 0, the ground state spm con- 
figuration is stable against a small perturbation. That 
is, the ground state phase will occur at least at very low 
temperatures. On the other hand, if 6* < 0, the ground 
state phase is absent at any non-zero temperature. 

To apply the above idea to our model, we must give 
special attention to the rotational axis for {Sl+ia} be- 
cause the reference system has a non- vanishing magneti- 
zation Mjf. For the following arguments, we separate 
each spin Si^i into parallel and perpendicular compo- 
nents: 

Sl = {Si,, ■ m)m 
S^^ — {Si^i X m) X m, 

where m = We consider two twisted sys- 

tems. One is a system in which {Sj^_^i j} are rotated 
around the axis that is parallel to the magnetization 
M'^ . We denote the minimum energy of this twisted sys- 
tem as Ej^{(j)). The other is a system in which {SL+i,i} 
are rotated around an axis that is perpendicular to ■ 
We also denote the minimum energy of this twisted sys- 
tem as Note that, in this twisted system, {5']'^} 
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FIG. 4: Excess energies ASf and 1\e\ for L x L x (L 4- 1) 
lattices for various spin concentrations: (a) x < xy and (b) 
X > xp. Open symbols represent ^Ej^ and filled symbols 
^e\. Symbols x in (a) represent the averages of those values. 

mainly change, but {^S^^ also change. Choices in the 
rotation axis are always possible in finite systems, even 
when X < xp because a non-vanishing magnetization 
(Af{f 7^ 0) exists in the Heisenberg model on a finite lat- 
tice. Of course the difference between E-^{(f)) and (c^) 
will diminish for L — > cx) in the range x < xp. The excess 
energies /S.Ej-^{4>) and Ai?|(0) in our model are given as 

Ai?f(0) = [Ei{cp)~E2l (3) 
A£;i(0) = [El{^)^Efl (4) 

with [• • • ] being the sample average. 

We calculated Ai?]^(0) and Ai?|((/)) for a common ro- 
tation angle of = 7r/2 in lattices of L < 14. Numbers of 
the samples are Ns ^ 1000 for L < 10 and Ns - 250 for 
L = 12 and 14. Hereafter we simply describe AEj^{7r/2) 

and Ai?|(7r/2) respectively as AE-^^ and AE^^. Figures 
4(a) and 4(b) respectively show lattice size dependences 
of AEj^ and Ai?| for x < xp and x > xp. We see that, 
for all x, AEj^ > AE-^^ and both increase with L. When 
X < xp, as expected, the difference between AEj^ and 

Ai?| diminishes as L increases. 

Now we discuss the stability of the spin configuration. 
First we consider the stability of {5']'^}, i.e., the stability 
of the FM phase. In the pure FM case {x ^ 1), S^.^ ^ 

and AEj^ gives the net excess energy for the twist of 
the magnetization. This is not the same in the case of 
S/t; 7^ 0. Because the twist in {S'j'j} accompanies the 
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FIG. 5: DifTerence in the excess energy AE[ = AEl - AEi FIG. 6: 
for L X L X (L + 1) lattice for various spin concentrations x. L x [L ■ 
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change in {S^^}, Ai?| does not give the net excess energy 
for the twist of {S'j'j}. For that reason, we consider the 
difference Ai?£ between the two excess energies: 

^El = ^El~^Ei. (5) 

If Ai?£ ^ cxo for L — > oo, the FM phase will be sta- 
ble against a small perturbation. We define the stiffness 
exponent 9^ of the FM phase as 

/^Elo^L'^\ (6) 

Figure 5 shows A£:£ for x > 0.80. We have 6*^ > for 
X > 0.85 and 0^ < for a; = 0.80. These facts show 
that, in fact, the FM phase is stable for x > xp ^ 0.84 
at r ~ 0. 

Next, we consider the stability of the transverse com- 
ponents {Sj^^}. Hereafter we call the phase with {Sj^^ ^ 
0} a SG phase. For x < xf, we may examine the stiff- 
ness exponent 6^'^ using either AEj^ or Ai?|. Here we 
estimate its value using an average value of them. For 
a; > xf, we examine it using AEj^. In this range of x, 
meticulous care should be given to a strong finite size 
effect ■ We infer that this finite size effect for x > xp 
is attributable to a gradual decrease in the magnetization 
M for finite L (see Fig. 2). That is, the magnitude of the 
transverse component \S^^\ will gradually increase with 
L, which will engender an additional increase of AE^ as 
L increases. This increase of \S[-^\ will cease for L oo. 
Consequently, we estimate the value of 6*^^ from the re- 
lations: 

(A£:|-|-A.Bf)/2 cx L^""^ for xKxp, (7) 
AEi/\S^\^ cx L^"*" for x > xp, (8) 

where \S-^\'^ = 1 - \M/xN\'^. Log-log plots of those 
quantities versus L are presented in Fig. 4(a) for x < xf 
and in Fig. 6 for a; > xf- We estimate 9^^ using data for 
L > 8 and present the results in the figures. Note that 
for X > 0.90, studies of bigger lattices will be necessary 
to obtain a reliable value of 9^"^ because AEj^ for L < 14 
is too small to examine the stiffness of {Sj^^}. 
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FIG. 7: Stiffness exponents O^'^ and 6^ for various spin con- 
centrations X. Here, we remove at a; = 0.95. 

Figure 7 shows stiffness exponents 9^ and 9^'^ as func- 
tions of X. As X increases, 9^'^ changes its sign from 
negative to positive at xsg = 0.175 ± 0.025. This 
value of Xsg is close to the percolation threshold of 
Xp ~ 0.137j23|. Above a::sGi ^^"^ takes almost the same 
value of 6*^*^ ^ 0.75 up to a; ^ 0.9. On the other hand, 
9^ changes its sign at a:F ^ 0.84 and increases toward 
61^ = 1 at a; = 1. A notable point is that 9^'^ > 
for a; > a:F- That is, a mixed (M) phase of the ferro- 
magnetism and the SG phase will occur for a; > a;F at 
T = 0. We could not estimate another threshold of x 
above which the purely FM phase is realized. 

IV. MONTE CARLO SIMULATION 

We next consider the magnetic phase at finite tem- 
peratures using the MC simulation technique. We make 
a MC simulation for x > 0.20. We treat lattices of 
LxLxL{L = 8 — 48) with periodic boundary con- 
ditions. Simulation is performed using a conventional 
heat-bath MC method. The system is cooled gradually 
from a high temperature (cooling simulation) . For larger 
lattices, 200000 MC steps (MCS) are allowed for relax- 
ation; data of successive 200000 MCS are used to calcu- 
late average values. We will show later that these MCS 
are sufficient for studying equilibrium properties of the 
model at a temperature range within which the RSG be- 
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FIG. 8: Specific lieats C in the 32 x 32 x 32 lattice for various 
spin concentrations x. 
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FIG. 9: Magnetizations M in the 32 x 32 x 32 lattice for 
various spin concentrations x. 



havior is found. Numbers Ng of samples with different 

spin distributions are A''^ = 1000 for L < 16, Ng = 500 for 
L = 24, Ns = 200 for L = 32, and = 80 for L = 48. 
We measure the temperature in units of Ji (fce = !)• 



A. Thermal and magnetic properties 

We calculate the specific heat C and magnetization M 
given by 



C 
M 



mis))]. 



(9) 
(10) 



Therein, E{s) and M(s)(= | J2i^i'^i\) represent the en- 
ergy and magnetization at the sth MC step, and A'' is the 
number of the lattice sites. Here (• • • ) represents an MC 
average. 

Figure 8 shows the specific heat C for various concen- 
trations X. For X > 0.90, C exhibits a sharp peak at a 
high temperature, revealing that a FM phase transition 
occurs at that temperature. As x decreases, the peak 
broadens. On the other hand, at a; ~ 0.85 a hump is ap- 
parent at a lower temperature; it grows with decreasing 
X. This fact implies that, for x < 0.85, another change 
in the spin structure occurs at a lower temperature. As 
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FIG. 10: Magnetizations M for a; = 0.80 in the L x L x L 
lattice. Open symbols indicate M in the cooling simulation 
and filled symbols indicate that in the heating simulation. 
Data at T = indicate those in the ground state given in 
Fig. 2. 



X decreases further, the broad peak at a higher temper- 
ature disappears and only a single broad peak is visible 
at a lower temperature. 

Figure 9 shows temperature dependencies of magneti- 
zation M for various x. For x = 1, as the temperature 
decreases, M increases rapidly below the temperature, 
revealing the occurrence of a FM phase. As x decreases, 
M exhibits an interesting phenomenon; in the range of 
0.78 < a; < 0.85, M once increases, reaches a maximum 
value, then decreases. We also perform a complemen- 
tary simulation to examine this behavior of M . That is, 
starting with a random spin configuration at a low tem- 
perature, the system is heated gradually (heating simu- 
lation) . Figure 10 shows temperature dependencies of M 
for = 0.80 in both cooling and heating simulations for 
various L. For T > O.lJi, data of the two simulations 
almost coincide mutually, even for large L. We thereby 
infer that M for T > O.lJi are of thermal equilibrium 
and the characteristic behavior of M found here is an 
inherent property of the model. For T < O.lJi, a great 
difference in M is apparent between the two simulations; 
estimation of the equilibrium value is difficult. We spec- 
ulate, however, that the lic^ating simulation gives a value 
of M that is similar to that in the equilibrium state be- 
cause the data in the heating simulation seem to concur 
with those obtained in the ground state. 

Figure 10 shows the remarkable lattice size dependence 
of M. For smaller L, as the temperature decreases, M 
decreases slightly at very low temperatures. The decrease 
is enhanced as L increases. Consequently, a strong size- 
dependence of M is indicated for T < O.lJi. These facts 
suggest that M for L ^ oo disappears at low temper- 
atures as well as at high temperatures. The next sec- 
tion presents an examination of this issue, calculating 
the Binder parameter. 
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B. Ferromagnetic phase transition 

The Binder parameter ijl at finite temperatures is de- 
fined as 



oiMf)!)Lwo 



(11) 



We calculate gL for various x. Figures ll(a)-ll(d) show 
g^s for X ^ xp^^- In fact, gL for x < xp exhibits a 
novel temperature dependence. As the temperature is 
decreased from a high temperature, gL increases rapidly, 
becomes maximum, then decreases. In particular, we 
see in Fig. 11(b) for x = 0.80 (Jl's for different L cross 
at two temperatures Tq and Tr (< Tq). The cross at 
Tc/Ji ~ 0.26 is a usual one that is found in the FM 
phase transition. That is, for T > Tq, gL for a larger size 
is smaller than that for a smaller size; for T < Tc, this 
size dependence in gL is reversed. On the other hand, the 
cross at Tr is strange: for T < Tr, for a larger size 
again becomes smaller than that for a smaller size. In- 
terestingly, the cross for different gL occur at almost the 
same temperature of Tr/Ji ~ 0.13. These facts reveal 
that, as the temperature is decreased to below Tr, the 
FM phase, which occurs below Tc, disappears. Similar 
properties are apparent for x = 0.79-0.82. 



C. Spin glass phase transition 

Is the SG phase realized at low temperatures? A con- 
vincing way of examining the SG phase transition is a 
finite size scaling analysis of the correlation length, ^l, 
of different sizes l[2^ [2^. Data for the dimensionless 
ratio £,l/L are expected to intersect at the SG tran- 
sition temperature of TsG- Here we consider the cor- 
relation length of the SG component of the spin, i.e., 
Si{= Si — m) with m as the ferromagnetic component 
of m = XiSij {xN). We perform a cooling simulation 
of a two-replica system with [Si] and {Ti} 27]. The SG 
order parameter, generalized to wave vector fc, q^^{k), is 
defined as 
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FIG. 11: Binder parameters ql for various x. The Tr for 
X — 0.82 was estimated by extrapolations of data obtained at 
higher temperatures. 



where = x, y, z. From this, the wave vector depen- 
dent SG susceptibility XSG(fc) is determinate as 



x^G{k)^xNYywm^)\ 



(13) 



The SG correlation length can then be calculated from 
1 . X5g(0) 



2sin(fc„iin/2) 

XSG(^min) 



1 



,1/2 



(14) 



where fcmin = (27r/L,0, 0). It is to be noted that, in the 
FM phase (m ^ for L — > oo), the FM component will 



interfere with the development of the correlation length 
of the SG component Si. Then in that case we consider 
the transverse components S^{= {Si x m) x m) in eq. 
(13) instead of Si. The correlation length obtained using 
S^ is denoted as 

We calculate ^l/^ or /L for 0.20 < a; < 0.90. The 
crosses for different L are found for 0.30 < a: < 0.90. Fig- 
ures 12(a)-12(c) show results of the temperature depen- 
dence of ^l/L for typical x. Assuming that the SG tran- 
sition occurs at the crossing temperature, we can scale 
all the data for each x (see insets). For x — 0.20, the 
crosses were not visible down to T j J\ — 0.02. However, 
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FIG. 12: The SG correlation length ^z, divided by L at dif- 
ferent X. Insets show typical examples of the scaling plot. 

we can scale all the data assuming a finite transition tem- 
perature of Tsq/Ji ^ 0.01. Thereby, we infer that the 
SG transition occurs for 0.20 x < 0.90. This finding 
is compatible with the argument in the previous section 
that > for 0.20 <x< 0.90. 

It is noteworthy that the SG phase transition for 
m 7^ is one in which the transverse spin components 
{S:^} order. Therefore we identify this phase transi- 
tion as a Gabay and Toulouse (GT) transition[2^ and 
the low temperature phase as a mixed (M) phase of the 
FM and a transverse SG. It is also noteworthy that, 
for X = 0.79 and x = 0.80, we estimate respectively 
Tsg/Ji = 0.10 ± 0.01 and Tgc/Ji = 0.098 ± 0.005, 
whereas respectively Tr/Ji = 0.15±0.01 and Tr/Ji — 
0.125 ± 0.0050. These facts suggest that, as the tem- 
perature is decreased, the SG transition occurs after the 
disappearance of the FM phase (TsG < 7r)- The differ- 
ence in transition temperatures of Tinv(= Tr) and Tsg 
were reported in Fco.tAIq.s 5]. However, further studies 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




PM 

- 


y - 


- 

Mixed 


1 

^ FM 


(FM+SG) 






'■)(, 







1 1 t t 

0.2 0.4 0.6 0.8 



I 



FIG. 13: The phase diagram of the dilute Heisenberg model. 
Four arrows indicate, from the left to the right, the percola- 
tion threshold Xp, the lower threshold of the SG phase xsg, 
the threshold of the ferromagnetic phase at finite tempera- 
tures xft, and the ferromagnetic threshold at T = 0, a;F. 



are necessary to resolve this point because the treated 
lattices of L < 20 for estimating Tsg are not sufficiently 
large. 



V. PHASE DIAGRAM 

Figure 13 shows the phase diagram of the model ob- 
tained in this study. It is shared by four phases: (i) the 
PM phase, (ii) the FM phase, (iii) the SG phase, and (iv) 
the M phase. A point that demands re-emphasis is that, 
just below the T = phase boundary between the SG 
phase and the M phase (xft < x < xp), the RSG tran- 
sition is found. This phase diagram is analogous with 
those observed in dilute ferromagnets Fej;Aui_a;0| and 
Eua;Sri_a;S0, 0. In particular, the occurrence of the 
mixed phase was reported in Fea;Aui_a;. 

We examine the low temperature spin structure. Fig- 
ures 14(a) and 14(b) represent the spin structure in the 
SG phase {x < xp). We can see that the system breaks up 
to yield ferromagnetic clusters. In particular, for a: < xp 
(Fig. 14(b)), the cluster size is remarkable. Therefore the 
SG phase for x < is characterized by ferromagnetic 
clusters with different spin directions. Figure 14(c) rep- 
resents the spin structure in the M phase {x > xp). We 
can see that a ferromagnetic spin correlation extends over 
the lattice. There are ferromagnetic clusters in places. 
The spin directions of those clusters tilt to different di- 
rections. That is, as noted in the previous section, the 
M phase is characterized by the coexistance of the fer- 
romagnetic long-range order and the ferromagnetic clus- 
ters with transverse spin component. The occurrence of 
ferromagnetic clusters at a; ~ xp are compatible with 
experimental observations [ll Q, H, H, ly] - 
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VI. CONCLUSION 

This study examined the phase diagram of a dilute 
ferromagnetic Heisenberg model with antiferromagnetic 
next-nearest-neighbor interactions. Results show that 
the model reproduces experimental phase diagrams of 
dilute ferromagnets. Moreover, the model was shown to 
exhibit reentrant spin glass (RSG) behavior, the most 
important issue. Other important issues remain unre- 
solved, especially in the RSG transition. Why does the 
magnetization, which grows at high temperatures, dimin- 
ish at low temperatures? Why does the spin glass phase 
transition take place after the disappearance of the ferro- 



magnetic phase? We intend the model presented herein 
as one means to solve those and other remaining prob- 
lems. 
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